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STORAGE  PROBLEMS  WHEN  DEMAND  IS  "ALL  OR  NOTHING" 

by 

D.  P.  Gaver  and  P.  A.  Jacobs 

Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA   9  3940 

0 .   ABSTRACT 

An  inventory  of  physical  goods  or  storage  space  (in 
a  communications  system  buffer,  for  instance)  often  experiences 
"all  or  nothing"  demand:   if  a  demand  of  random  size   D   can  be 
immediately  and  entirely  filled  from  stock  it  is  satisfied, 
but  otherwise  it  vanishes.   Probabilistic  properties  of  the 
resulting  inventory  level  are  discussed  analytically,  both 
for  the  single  buffer  and  for  multiple  buffer  problems . 
Numerical  results  are  presented. 


1 .   INTRODUCTION 

The  usual  storage  or  inventory  problems  involve  demands 
imagined  to  occur  randomly,  and  to  be  capable  of  reducing  any 
available  stock  to  zero,  or  even  beyond  when  backordering 
is  permitted.   Yet  in  many  situations  at  least  one  component 
of  total  demand  is  "all  or  nothing; "  that  is,  it  reduces 
inventory  only  if  it  can  be  entirely  satisfied  by  the  inventory 
present,  and  otherwise  seeks  another  supplier.   Here  are 
examples . 
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(a)   A  manufacturer's  warehouse  is  filled  with  a  certain  item 
at  the  beginning  of  the  selling  season;  let   I   denote 
the  initial  inventory.   Suppose  that  demands  occur  as 
follows:   a  message  is  sent  requesting  that   D.   items 
be  shipped  from  inventory,  but  only  if  the  entire  order 
can  be  filled.   That  is,  the  demand  is  satisfied  if   D,  <_  I , 
in  which  case  inventory  level  is  reduced  to   1(1)  =  I  -  D.  ; 
while  if   D,  >  I   the  inventory  remains  unchanged 
and   1(1)  =  I.   Allowing  for  no  replenishment,  the  second 
demand,  of  size   D~ ,  interacts  with  inventory   1(1), 
so  that  it  is  filled  if   D2  <_   1(1),  but  is  not  placed  if 
D2  >  1(1).   The  process  continues  along  these  lines  until 
the  selling  season  is  over  and  there  are  no  more  demands. 

(b)  A  buffer  storage  device  used  to  contain  messages  prior  to 
their  batch  transmission  has  capacity   I.   Messages  of 
length   {D^,  i  =  1,2,...}   approach  the  buffer  successively, 
and  are  admitted  on  an  "all  or  nothing"  basis,  just  as 

was  true  of  demands  for  physical  inventory  in  (a)  above. 
Once  again  rejection  will  occur,  and  more  frequently  to 
large  demands  (messages)  than  to  short  ones. 

(c)  A  system  of  many  buffer  storage  devices  is  used  to  contain 
messages  prior  to  their  batch  transmission.   Each  buffer 
has  capacity   I.   Messages  of  length   {D. ,  i  =  1,2,...} 
approach  the  device  and  are  successively  admitted  to  the 
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first  buffer  until  there  is  a  demand  that  exceeds  its 
remaining  capacity.   The  first  buffer  is  left  forever 
and  the  demand  that  exceeds  the  first  buffer,  plus 
successive  demands,  apply  to  the  second  buffer  until  one 
occurs  that  exceeds  the  remaining  capacity.   This  demand 
then  applies  to  the  third  buffer,  and  so  on.   As  a 
result  there  will  be  some  unused  capacity  in  each 
buffer. 

In  Section  2  we  will  discuss  some  models  for  the 
situations  in  examples  (a)  and  (b) .   We  will  compute  such 
quantities  as  the  distribution  of  the  amount  of  inventory 
left  at  some  time  t  and  the  distribution  of  the  times  of 
successive  unsatisfied  demands. 

In  Section  3  we  will  consider  a  model  for  example  c. 
We  will  derive  equations  for  the  limiting  distribution  of 
used  capacity  of  a  buffer  and  the  expected  used  capacity  of 
a  buffer.   It  seems  to  be  difficult  to  obtain  simple  analytic 
solutions  to  these  equations,  but  we  will  present  certain 
illustrative  numerical  results. 


2.   THE  ONE-BUFFER  INVENTORY  PROBLEM 

Suppose  that  demands  for  available  stock  occur 
according  to  a  compound  Poisson  process:   if  N    is  the 
number  of  demands  that  occur  in   (0,t],  then   'tNt/*'t  2.   0} 
is  a  stationary  Poisson  process  with  rate   A;  the  sizes  of 
successive  demands   {D.}   are  independent  with  common 
distribution   F.   Assume  that  there  are  no  replishments 
of  inventory.   Let   {lt/*  t  >^  0}   denote  the  stochastic 
process  describing  available  inventory  at  time   t,  and 
let   {I (n) ;  n  =  0,1,...}   be  the  stochastic  process  of 
available  inventory  following  the  nth  demand.   It  is 
apparent  from  our  assumptions  that  both  {*±.}      and   (I(n)} 
are  Markov  processes. 

2.1.   Functional  Equations  for  the  Amount  of  Available  Inventory 
Let 


-si. 

(2-D  <Ks,t)  =  E[e    Z] 


be  the  Laplace  transform  of  the  available  inventory  at  time  t 
Similarly,  let 

iMs,n)  =  E["sl(n)]  . 

Properties  of  the  available  inventory  can  be  studied  in  terms 
of   0   and  ty;    we  begin  by  deriving  an  equation  for   <J>. 

4 


Observe  that  if  one  conditions  on   I  ,  then  in   (t,  t +  dt) 
either  no  change  in  inventory  occurs,  an  event  of  probability 
1  -  AF(I.)dt  +  o(dt),  or  a  depletion  of  amount   x   occurs 
with  probability   A  dt  F(dx)  ,  x  <_  I  .   Thus 


(2.2)   E[exp(-slfc+dt) |lt] 

-si  ^t 

=  e    t[l  -  XF(I.)dt]  +  J      exp[-s(It  -  x)  ]  Adt  F(dx)  +  o(dt) 

0 


Now   take   expectations  with   respect   to      I 


(2.3)       4>(s,    t  +  dt) 


-si.  *t 


=    <J)(s,t)  -  A{E[e         fcF(I    )]-E[/      exp{-s(I    -x)}  F(dx)  ]}  dt  +  o(dt) 

z  0  z 


After  subtraction  of  <J>(s,t)   from  both  sides,  division  by   dt, 
and  allowing   dt  ■>  0   we  find  that   cf)   must  satisfy  the 
equation 


I 
0 
An  analogous  argument  shows  that 


-si    Xt 
(2.4)  M  =  AE  j"e    t    J         (eSX-l)  F(dx) 


,  v   l(n) 
(2.5)      ip(s,  n+1)  =  iMs,n)  +  E  [e  sl(n)    /     (esx  -  1)  F(dx)]. 


Differentiation  with  respect  to   s   at   s  =  0 ,  or  a  direct 
conditional  probability  argument,  now  produce  equations  for 
E[I  ]   and   E[I(n)  ]  : 

d  Xt 

(2.6)  -£r  E[I  ]  =  -\E[J      x  F(dx)  ] 

r         0 


and  I(n) 

E[I(n+l)]  =  E[I(n)]  -  E[/     x  F(dx)]  . 

0 


In  general  no  explicit  solutions  for  the  expected 
values  are  available,  but  a  simple  upper  bound  results  from 
rewriting  (2.6)  as  follows. 

d  't 

(2.7)  -S-  E[I  ]  =  -XE[I.  /   JL  F(dx)] 

at     t  r  0   Xt 

>  -AE[lt  F(It)] 

>  -AF(I)  E[It]  , 

from  which  one  sees  that 


(2.8)  E[It]  >  lexp[-AF(I)t] 

and  similarly 

E[I(n)  ]  >  l[l  -  F(I)  ]n  , 

so  the  expected  available  inventory  declines  by  at  most  an 
exponential  rate. 


2.2.   Explicit  Solution  When  the  Demand  Distribution  is  Uniform 

Although  Equation  (2.4)  seems  to  be  quite  intractable 
for  most  demand  distributions , it  can  be  solved  completely 
when   F   is  uniform: 


x 

c 


F(x)  = 


0  <  x  <  c, 


C  >  X 


and   c  >  I.   In  this  case  (2.4)  can  be  expressed  as 


(2.9) 


H   =  «  [*"Slt  /^ 


=  AE 


0 
-si 


,  sx    , .  dx 
(e   "  1)  — 


-si 


] 


1  -  e 


sc 


_  A  ["1  -  4)1    A 
c  |_   s      c 


8s 


In  other  words  $      satisfies  a  first-order  (quasi) 

linear  partial  differential  equation  with  initial  condition 

-si 
<t>(s,0)  =  e    .   The  usual  procedure  for  solution,  Sneddon 

(1957) /  requires  solution  of  two  ordinary  differential 

equations  selected  from  among 


(2.10) 


dt      -ds     dj> 

1      (A/c)    (A/c)  [(l-<J>)/s] 


we  find  from  the  first  and  last  two  that 


A  1  -  £ 

(2.11)        s  +  -  t  =  C]_  ,        — —  -  c2  , 


so  a  general  solution  is  given  by 

(2.12)  g(s  +  |t,  *-=-*)-  0  - 

g  being  a  function  to  be  determined.   Now  the  specified 
initial  condition  stipulated  that  at   t  =  0 


(2.13)  *  -  e"SI  =  0  =  g(s,  L^_* 


SO 


(2.14)  i-^£  -  *  "  fSl  =  0 


which  specifies  <$>      at   t  =  0.   But  for   t  positive  we 
replace   s,  the  first  argument  of  g   at   t  =  0   by 
s  +  (A/c)t   to  obtain  the  solution 


(2  15)  1  -  4>(s,t)  _  1  -  exp[-(s  +  (X/c)  t)I]  =  Q 

s  s  +  (X/c) t  ' 


which  gives  the  desired  transform.   Passage  to  the  limit  as 
s  ->■  0   in  (2.15)  shows  that 


(2.i6)  s'v  ■  -  "(x/i';(At/c)I1  • 


This  formula  can  also  be  derived  by  first  finding  an  expression 
for  the  kth  moment  of   I  ,  and  then  employing  a  Taylor  series 
argument. 

In   order    to    invert    the    transform   in    (2.15)    note    that 


(2.17)         /   e"sx  P(I,    >   x>dx  =   1   ~    ^(s-t}    =   1  ,-  exp[-(s+  Ut/c))I] 

p.  u  S  S     t      lA/C)  t 

which  is  the  transform  of  a  truncated  exponential  distribution. 
Thus  by  the  unicity  theorem  for  Laplace  transforms 


(2.18)     P{I   >  x}   = 


exp[-(At/c)x]       0  <  x  <  I, 


I  <  x 


Note  that  the  distribution  of   I    is  absolutely  continuous 
in  the  interval   (0,1)   but  that  there  is  a  jump  at   I 
corresponding  to  the  occurrence  of  no  demand  less  than  or 
equal  to   I   in   (0,t]: 

(2.19)  P{ifc  =  1}  =  exp[-  Xt(I/c)]  . 


2.3.   The  Expected  Number  of  Satisfied  Demands 

Supposing  that  an  initial  inventory,  or  storage  capacity, 
I   prevails,  it  is  of  interest  to  compute  the  probability  that 
a  demand  is  satisfied,  and  the  expected  number  of  demands 
satisfied  in  an  interval  of  length   t.   First  notice  that  if  a 
demand  of  size   D(t)   appears  at  time   t,  at  which  moment   I 
is  available,  then 

P(D(t)  <  Itllt>  =  F(It) 

is  the  conditional  probability  that  the  demand  is  satisfied. 
When  F  is  uniform,  as  is  presently  true,  we  may  remove  the 
condition  to  find  that 

P(D(t)  <  I)  =  E[F(IJ]  =  E  f  -^  1  =  1  -  exp[-(At/c)I] 
-  t         t       |_  c  J         At 

If   S(t)   is  the  number  of  demands  satisfied  during  the  time 
interval   (0,t],  then  since  demands  arrive  according  to  a 
Poisson  process  with  rate   A, 


t  t 

(2.20)    E[S(t)]  =  A  /   E[F(I  )  ]  du  =  A  /   1  -  exp [- (Au/c) I ]  du 

0        u  0  Au 

=  Y  +  £n(A|i)  +  E  (Atl} 
c      i   c 

where  E1  ( • )   is  an  exponential  integral;  Abramowitz  and 
Stegun  (1965),  and   y  =  0.5112...   is  Euler's  constant. 
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2.4.   The  Time  of  the  First  Unsatisfied  Demand  and  the  Amount 
of  Unused  Inventory  at  That  Time 

As  before   F   is  the  common  distribution  function  of 
the  successive  demands.   Now  let   t   be  the  time  of  the  first 
unsatisfied  demand.   Then 


P{  T  >  t|N.  =  n}  =  P(Dn  <  I,  D0  <  I  -  Dn  ,  .  .  .  ,  D  <  I  -  D- D   J 

t  1—2—     1 '       n—     1        n-1 

=  F(n)(I) 


where   F     denotes  the  nth  convolution  of   F  with  itself 
Hence 


n 


(2.21) 


P(x  >  t>  =  I      e"At  i*£L.  F<n>( 
n=0 


n: 


Explicit  expressions  for  the  distribution  of  t   can 
be  obtained  in  some  cases.   If   F   is  uniform  on   [0,c]   with 
c  >  I,  then 


(2.22) 


P{t  >  t}  =  e"At  I, 


Atl\ 


1/2 


where   Iq(z)   is  a  modified  Bessel  function  of  the  first  kind 
of  the  zeroth  order.   In  this  case 


(2.23)         e[t]  =  y  exp{l/c}  =  T  exp{l/2E[D]} 
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If  F  is  exponential  with  mean   1/y,  then 


n=0  k=n 


and 


(2.25)         E[x]  =  I  [1  +  ul]  =  |  [1  +  g-j^y] 


Note  that  if   I   is  small  relative  to  E[D],  then  the 
expected  time  to  first  unsatisfied  demand  when   F   is  exponential 
will  be  greater  than  the  expected  time  when   F   is  uniform.   How- 
ever for   I   large  relative  to  E [D]  the  expected  time  for   F 
exponential  will  be  less  than  the  expected  time  when   F   is 
uniform. 

Let  Y   be  the  amount  of  inventory  present  at  the  time 
of  the  ntn  unsatisfied  demand.   Then  for   0  <  a  <  I 


(2.26)  P{Y   >  I  -  a)  =  /   R(dy)  F(I  -  y) 

0 
where 

oo 

(2.27)  R(y,  =  I      F(n) (y) 

n=0 
and 

(2.28)  F(I  -  y)  =  1  -  F(I  -  y) . 

Again  explicit  expressions  for  the  distribution  of  Y 
can  be  obtained  for  some  distributions  F.  If  F  is  uniform 
on   [0,c]   for   c  >  I,  then 
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(2.29)         P(Y1    >    I    -   a)    =    1    -    i1    ~   a)    exp{|  a}     . 


If      F      is   a   truncated  exponential 


(2.30)         f(x)    = 


,  -yx 

1    -   e    H 

1    -  e    H 


x  <    I    , 


x   >   I    , 


then 

(2.31)  P{Yl  >  I-a}  =  1-  [e~ya-e~yI]  [l-e"yI]_1  exp{pa  [l-e-^1]'1} 

If  F  has  an  exponential  distribution  with  mean   1/y ,  then 

(2.32)  P{Yl  >  I  -  a}  =  e"y(I"a)  . 

In  this  last  case  the  distribution  function  of  Y   can  be 

n 

computed  by  induction  quite  easily  and 

(2.33)  P{Y   >  I  -  a}  =  e"ny(I_a)  . 

n  — 

Hence  when      F      is   exponential 


(2.34)  E[Y    ]    =  -i-    [1    -  e'nyI] 

n  ny 


In  principle  similar  results  can  be  obtained  for  other  distri- 
butions, but  we  have  found  no  simple  expressions. 
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2.5.   Inventory  Costs  and  Policies 

There  are  at  least  three  monetary  quantities  which  affect 
the  profitability  of  an  inventory  policy  over  a  fixed  interval  of 
time   (0,t]:   the  selling  price,  p;  the  storage  cost,  a;  and 
the  cost  of  lost  demands,  b.   If  the  storage  cost   a   is 
charged  just  on  the  basis  of   I   (something  like  warehouse  size) 
then  the  total  expected  profit  in   (0,t]   is 

Z(I)  =  p(i  -  E[It])  -  al  -  bS(t) 

-1 


=  (p  -  a)  I  -  p 


(|  t)    [1  "  exp[-(At/c)I]] 


-  b{Y  +  m(±±   I)  +  E,  (^   I)} 


for  the  case  of  uniformly  distributed  demands;  see  ((2.16)  and 
(2.20)).  One  can  numerically  find  the  maximum  expected  profit 
for  this  case;  nothing  explicit  seems  to  be  available. 


14 


3.   THE  MANY-BUFFER  STORAGE  PROBLEM 

In  this  section  we  will  study  a  model  for  the  situation 
of  example  (c)  in  section  1.   Messages  are  successively  admitted 
to  the  nth  buffer  until  there  is  a  message  length  that  exceeds  the 
remaining  capacity  of  the  buffer.   The  total  amount  of  this  message 
is  put  in  the  (n+l)st  buffer  and  the  nth  buffer  is  left  forever. 
Successive  messages  are  then  put  in  the  (n+l)st  buffer  until  there 
is  a  message  whose  length  exceeds  the  remaining  capacity  of  the 
(n+l)st  buffer;  this  message  is  put  in  the  (n+2)nd  buffer  and 
so  on . 

Let   I   denote  the  common  capacity  of  the  buffers  and 

D.   denote  the  length  of  message  i.   Assume   {D.}   is  a  sequence 

of  independent  identically  distributed  random  variables  with 

distribution   F  having  a  density  function   f   such  that 

f(x)  >  d  >  0   for   x  6  [0,1].   Let  R(x)  =  Y"  n  F(n) (x)   be 
—  ^n=U 

the  renewal  function  associated  with   F.   If  F(I)  <  1,  then 
we  will  assume  that  an  incoming  message  to  the  currently  used 
nth  buffer  of  length  greater  than   I   is  sent  to  the  (n+l)st 
buffer;  when  it  cannot  fit  into  the  (n+l)st  buffer,  then  it  is 
"banished,"  i.e.  sent  to  some  other  set  of  buffers.  The  next 
message  however  will  try  to  enter  the  (n+l)st  buffer.   If  this 
message  has  length  greater  than   I   it  is  banished  and  the 
following  message  will  try  to  enter  the  (n+1) st  buffer;  all 
messages  of  length  exceeding   I   will  be  banished  until  one 
appears  that  is  smaller  than   I   and  it  will  be  the  first 
entry  in  buffer  (n+1) . 
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This  model  has  been  studied  for  demand  distributions   F 

with   F(I)  =  1   by  Coffman  et  al.  (19  78) .   Their  approach  was 

to  study  the  Markov  process  describing  the  total  amount  of 

inventory  or  space  consumed  in  successive  buffers  or  bins. 

Here  we  study  the  process   {L  },  where   L    is  the  size  of  the 

n  n 

demand  that  first  exceeds  the  remaining  capacity  of  the  nth 
buffer;   {L  ;  n  =  1, 2 , . . . }   is  a  Markov  process.   Let 


K(x,  [0,y])  =  P(Ln+1  £  ylLn  =  *} 


Note  that 


P{L1  £  y}  =  K(0,  [0,y]) 


is  the  same  as  the  sum  of  the  forward  and  backward  recurrence 
times  at  time   I   for  a  temporal  renewal  process  with  inter- 
renewal  distribution  F.   Thus  for  y  <  I 


(3.1)      H1(y)  =   P{L1  £  y>  =  /    R(dz)  [F(y)  -  F(I-z)] 

I-y 


Note  that  for  y  <  I 


(3.2)   K(x, [0/Y])  = 


I-x 


I-x-y 


R(dz)  [F(y)-F(I-x-z)  ] 


if   x  <  I-y; 


I-x 

/    R(dz) [F(y)-F(I-x-z) ] 

0 

if   I-y  <  x  <  I; 


/    R(dz)  [F(y)-F(I-z)  ] 
I-y 

if   x  >  I 
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Hence 


[R(I-x)  -  R(I-x-y)]  F(dy)      if   x  <  I-y, 


R(I-x)  F(dy) 


if   I  >  x  >  I-y, 


(3.3)   K(x,dy)  = 


|R(y)  F(dy)  -  /   R(dz)  f(y-z)  +  R(dy)F(y) 
0 


[R(D  -  R(I-y) ]  F(dy) 


if   x  =  I-y  , 
if   x  >  I  . 


Note  that  for  some   0  <  a  <  b  <  I,  there  exists  a   6  >  0   such 
that  for  all   x 

2 
K  (x,dy)  2  5      for     y  €  [a,b] 

where   K  (x,dy)  =  /~  K(x,dz)  K(Z/dy).   Hence  hypothesis  D' 
on  page  197  of  Doob  (1952)  is  satisfied.   Thus,  if 


Kn(x,A)  =  P{L1+n  i    All^  =  x> 


for  all  Borel  subsets   A,  then 


(3.4) 


,n 


lim  K  (x,A)  =  H(A) 

n  ■*  °° 


exists  and  further  the  convergence  is  geometric 

|Kn(x,-A)  -  H(A)  I  <  ay" 

for  some  positive  constants   a   and   y,   y  <  1   for  all   A, 
Now  let 
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H  (x)  =  P{L   €  [0,x] |Lft  =  0} . 
n         n  u 


Then  a  renewal  argument  can  be  used  to  show  that  for   x  £  I 

I 
(3.5)       Hn+i(x)  =  /   Hn*R(dy)  [F(x)  -F(I-y)] 

I-x 


+  [1-H  (I)]  /    R(dy)  [F(x)  -F(I-y)] 
I-x 


Taking  limits  as   n  -*  °°   it  is  seen  that  the  distribution 
H(x)   satisfies  the  following  equation  for   x  <  I: 


I 
(3.6)       H(x)  =  /    H*R(dy)  [F(x)  -  F(I-y)] 

I-x 

I 
+  [1-H(D]  /    R(dy)  [F(x)  -  F(I-y)]  . 
I-x 


Equations  (3.1)  and  (3.6)  can  be  simplified  for  certain 
specific  distributions   F. 

A .   Exponential  Demands . 

For  the  exponential  distribution  with  mean  1  and 
x  <  I   the  equations  are 


(3.7)  H1(x)  =l-ex-xeX 


and 

-x  -x  ,x 

(3.8)       H(x)  =  xe    H(I)  +  H1(x)  -  e  x  /   H(I-x+u)du 
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B.       Uniform   Demands 

For    the    uniform  distribution   on       [0,c]      with      c    >    I 
they   simplify    to 

(3.9)  H-^x)    =   exp[|    (I-x)]    -    (1    -  |)    exp(|  I) 
and 

I-x 

(3.10)  H(x)    =   i  exp[i    (I-x)]    /         exp(-   ^  u)    H(u)du 

c  c  0  c 

+   *-  H(I)    -   i  exp(^   I)     (1    -   *-)    /      exp(-   ±  u)    H(u)du 


+    [1    -   H(I)  ]    H-^x)     , 


for   x  <_  I .   Similar  expressions  hold  for   x  >  I,  but  they  are 
unimportant  in  the  present  context. 

Equations  (3.6),  (3.8)  and  (3.10)  do  not  seem  to  yield 
explicit  answers.   As  a  result  we  have  solved  (3.8)  and  (3.10) 
numerically  by  iteration  using  the  system  of  equations 


x 

(3.11)   H   , (x)  =  xe  x  H  (I)  +  Hn  (x)  -  e  x  /   H  (I-x+u)du 
n+l  n       1  jL    n 


with   H,   as  in  (3.7)  and 


l      l        I_x       l 
(3.12)   H  ,1 (x)  =  i  exp[-  (I-x)]  /    exp(-±u)  H  (u) 
n+l      c  c  i  c     n 


+  -  H  (I)  -  -  exp(-  I)  (1  -  -)  /   exp(-  -  u)  H  (u)du 
en       cc        ci^c     n 


+  [1  -  H  (I)  ]  H.  (x) 
n      1 


with   H,   as  in  (3.9).   For  the  cases  carried  out  the  convergence 
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is  rapid;  after   n  =  5   iterations  very  little  change  is  noted 
and  convergence  has  occurred,  for  most  practical  purposes. 

Next  let   Y    be  the  amount  of  storage  space  used  in 
the  nth  bin;  the  distribution  of   Y    is  denoted  by   G  (x) ,  and 


G(x)  =   lim   P(Y   <  x)  =   lim   G  (x) 

n  —  n 

n  ->  °°  n  •*•  °° 


is  the  long-run  distribution.   By  probabilistic  arguments 
and  (3.4) 

x  * 

(3  13)   G(x)  =  /   H*R(dy)  F(I-y)  +  [l-H(I)]  /   R(dy)  F(I-y) 
0  0 

where   F(I-y)  =  1  -  F(I-y)   and  the  long  run  average  expected 
capacity  of  a  bin  that  is  actually  used  is 

I 
A  =  /   x  G(dx)  . 
0 

For  the  case  in  which   F   is  exponential  with  unit  mean 


(3.14)      A  =  I  -  [l-H(I)][l-e  I]  -  e  Z  /   eX  H(x)dx  . 

0 


For  the  case  in  which   F   is  uniform  on   [0,c]   with   c  >_  I 

I  1      I       1 

(3.15)   A  =  -2  /   H(u)d   +  exp(-  I)  /   exp(-  -  u)  H(u)du 

0  c     0         C 

+  H(I)  [21  -cexp(-  I)  +  c]  +  [-1  +  c  exp(-  I)  -  c] 
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Numerical  solutions  were  obtained  for  equations  (3.14)  and 

(3.15)  by  first  computing  the  probabilities    H  (x) ,  n  =  1,2,..., 10 

iteratively  from  (3.7)  and  (3.11)  for  the  exponential  demand 

case,  and  from  (3.9)  and  (3.12)  for  the  case  of  uniform  demands. 

Our  technique  was  simply  to  discretize  x:  x.  =  jh,  h  =  I/N, 

N   beinq  the  number  of  x-values  at  which   H  (x)   is  evaluated 

n 

(values  of   N   from  200-1200  were  utilized  in  order  to  obtain 

two-significant  digit  accuracy) .   The  integrals  were  then 

approximated  by  a  summation,  i.e.  Simpson's  rule.   Having  the 

values  of   H  (x.)  it  is  possible  to  calculate  those  of 

H  ,, (x.) ,  and  from  these  the  values  of   G  (x)   and  the  mean 
n+1   j  n 

usage,  E [Y  ],  may  be  calculated  by  numerical  integration.   In 
the  case  of  exponential  demand  very  simple  upper  and  lower 
bounds  were  obtainable;  such  bounds  were  not  tight  enough  to 
be  useful  for  the  uniform  case. 

The  following  table  summarizes  the  numerical  results. 
We  have  compared  demand  distributions  that  result,  as  nearly 
as  possible,  in  the  same  probability  that  an  initial  demand 
on  an  empty  bin  will  be  rejected.   We  have  tabulated  the 
expected  level  to  which  the  bin  is  filled.   It  is  interesting 
that  the  limited  bin  occupancy  is  0.75  when  a  uniform  demand 
over  the  range  of  the  bin  size  is  experienced.  This  result  has 
been  obtained  analytically  by  Coffman  et  al .  (1978);  in  that 
paper  simple  and  elegant  analytical  expressions  for   G   and   H 
also  appear  for  this  case.   The  considerable  similarity  of 
the  numbers  in  the  rows  of  the  table  is  notable;  apparently  the 
long-run  bin  occupancy  is  only  slightly  larger  than  is  that 
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of  the  first  bin,  and  the  occupancy  experienced  for  uniform 
demand  is  only  slightly  larger  than  for  exponential.   Further 
investigations  to  examine  the  reasons  for  this  insensitivi ty 
would  seem  to  be  of  interest. 

Expected  Fraction  of  Bin  Filled 
<fn  =  E[Yn]  .  I) 

Rejection  Probability     Exponential  Demand    Uniform  Demand 


F(I) 
0.00 
0.05 
0.10 
0.15 
0.20 
0.25 


fl 

f 

oo 

fi 

f 

oo 

- 

- 

0.76 

0.75 

0.74 

0.75 

0.74 

0.74 

0.69 

0.70 

0.72 

0.72 

0.65 

0.66 

0.68 

0.69 

0.60 

0.62 

0.64 

0.66 

0.56 

0.58 

0.60 

0.62 
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